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I. INTRODUCTION 

Excitons (electron-hole bound states) within quantum 
dots (QD's) have already attracted much interest in the 
field of quantum computation and have formed the basis 
of several proposals for quantum logic gatesi The en- 
ergy shift due to the exciton-exciton dipole interaction 
between two QD's gives rise to diagonal terms in the in- 
teraction Hamiltonian, and hence it has been proposed 
that quantum logic may be performed via ultra-fast laser 
pulses^ However, excitons within adjacent QD's are also 
able to interact through their resonant (Forster) energy 
transfer, 3 some evidence for which has been obtained ex- 
perimentally in a range of systems4ui*§*2i As is shown 
below, this resonant transfer of energy gives rise to off- 
diagonal Hamiltonian matrix elements and therefore to a 
naturally entangling quantum evolution. It is proposed 
here that this interaction may be observed in a straight- 
forward way through the observation of anticrossings in- 
duced in the coupled dot energy spectra by the applica- 
tion of an external static electric field. We also show that 
the off-diagonal matrix elements can be made sufficiently 
large to be of interest for excitonic quantum computa- 
tion^ 

The first studies of Forster resonant energy transfer 
were performed in the context of the sensitized lumines- 
cence of solids Here, an excited sensitizer atom can 
transfer its excitation to a neighbouring acceptor atom, 
via an intermediate virtual photon. This same mecha- 
nism has also been shown to be responsible for exciton 
transfer between QD's^ and within molecular systems 6 
and biosystems^ (though incoherently, as a mechanism 
for photosynthesis), all of which may be treated in a sim- 
ilar formulation. Thus, the results reported here can be 
expected to apply not only to QD systems, but also to 
a wide range of nanostructures where Forster processes 
are of primary significance. 

In this paper, we consider two coupled generic QD's, 
each modelled by a simple potential which is given by 
infinite parabolic wells in all three (x, y, z) dimensions. 
This potential profile allows an analytical expression for 



the interdot Forster coupling in the dipole-dipole approx- 
imation to be derived and, although we expect it only to 
allow for qualitative predictions about real observations, 
similar models have been successfully used in the liter- 
aTiUre I ii 4 i2 i i3 Excitations of each dot are assumed to be 
produced optically, and we neglect tunneling effects be- 
tween the two coupled dots (see the Appendix for more 
details). We shall restrict ourselves to small, strongly 
confined dots, with electron and hole confinement ener- 
gies ~ 100 meV and low temperatures (T < 5 K; see 
Sec. llVfl . We shall therefore consider only the ground 
state (no exciton) and first excited state (one exciton) 
within our model, as they will be energetically well sep- 
arated from higher excitations. These two states define 
our qubit basis as |0) and |1), respectively. 

Various techniques exist for growing QD's in the lab- 
oratory, of which the Stranski-Krastanow methodM^ is 
possibly the most promising for the realization of a con- 
trollably coupled many-dot system. In this growth mode 
a semiconductor is grown on a substrate which is made of 
a different semiconductor, leading to a lattice mismatch 
between the layers. Under certain growth conditions, 
dots form spontaneously due to the competing energy 
considerations of dot surface area, strain, and volume. If 
a spacer layer of material is then grown above the first 
dot layer and then a second dot layer deposited, a ver- 
tically correlated arrangement of dots can be made^ 
Two such stacked dots could then form our interacting 
two-qubit system with materials, growth conditions, and 
spacer layer size tailored to give suitable electronic prop- 
erties and inter-dot Coulomb interactions. 



II. QUANTUM DOT MODEL 

A. Single-particle states 

Many varying approaches to the calculation of elec- 
tron and hole states i n QD's have been put forward in 
Refs. l2lllll2ll3ll7ll8llT9ll20ll2ll the choice of which de- 
pends on the final aim of the work. The aim here is 
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to provide a clear and simple illustration of how to ex- 
perimentally observe resonant energy transfer between a 
pair of QD's and how to exploit this inter-dot interac- 
tion to perform quantum logic. Therefore, we shall con- 
sider oneofthe most basic models, similar to those in 
Refs. lllll2ll.il which treats the conduction- and valence- 
band ground states as those of a three-dimensional in- 
finite parabolic well. This simple model can provide us 
with analytical expressions for the energies and wave- 
functions of the single-particle states and for the dipole- 
dipole interaction between two dots. We do not expect 
that this model will predict the precise values of experi- 
mental measurements; however, it is expected to indicate 
qualitative trends in real observations. 

In the effective mass and envelope function approxi- 
mations22i2£ the Schrodinger equation for single particles 
may be written as 



EgbAs = 1.425 eV 




ffi(r)&(r) 



h 2 f 1 
2 Vm* 



V + Vi(r) 



where i — e,h for electron or hole, Vi(r) is the dot con- 
finement potential which accounts for the difference in 
band gaps across the heterostructure, and m* is the ef- 
fective mass of particle i. Here, <fii(r) is the envelope 
function part of the total wavefunction: 



FIG. 1: The parameters Cj are chosen from electron and hole 
confinement potentials V e and Vh by taking (l/2)cjr 2 — V%, 
</>i(r) = Ei<j)i{r), for j — x,y,z, and i — e,h. This matches a well of depth V5 
with the parabolic potential (l/2)cjj 2 at a width r from the 
(1) dot centre. All parameters are chosen to be consistent with 

V e + V h + E qap = 1.425 eV, the GaAs 



&(r) = ^(r)^(r). 



(2) 



The envelope function describes the slowly varying con- 
tribution to the change in wavefunction amplitude over 
the dot region, and the physical properties of the single- 
particle states can be derived purely from this contri- 
bution. Ui(r) is called the Bloch function and has the 
periodicity of the atomic lattice. Its consideration is vi- 
tal when describing the interactions between two or more 
particles. 

In the approximate analytical model, a separable po- 
tential comprising infinite parabolic wells in all three di- 
mensions represents the QDi£i 
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V(x,y, z) = -Ci. x x + -Ci. y y + -c^ z z 



(3) 



GaAs spacer layers, i.e. 
band-gap. 
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in the x- direction with analogous expressions for y and z. 
The integer n — (0, 1,2,3, ...) labels the quantum state, 
with energy E n = (n + l/2)hw x , the H n 's are Hermite 

polynomials, and d x = {hj ' \/m*c x ) 1 = [h/ {m*uj x )) 1 ^ 2 . 
We are interested only in the ground state solutions of 
each well, so our envelope function is given by: 
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where the frequency u>ij = y/a.j/m* for j = x,y,z 
(see Fig. Q). Hence, the Schrodinger equation (Eq. 
is also separable and provides simple product solutions 
for the electron and hole states. The envelope functions 
are therefore given by 



with energy Eq = ^h(u> x + ui y + uj z ). The choice of con- 
stants Cj and hence dj will be different for changing con- 
finement potentials and particle masses, and so will de- 
pend upon the energies of the system under consideration 
and whether the particle is an electron or hole (see FigOJ. 



(4) 



for the parabolic confinement (r = (x, y, z)). We now 
drop the subscript i but remember that due to their dif- 
fering effective masses electrons and holes may take dif- 
ferent values for the parameters defined throughout this 
paper. 

The solutions to the one-dimensional Schrodinger 
equation for the potential form of Eq. [21 are given by: 



B. Excitons and Coulomb integrals 

The excitation of an electron from a valence band state 
to a conduction band state leaves a hole in the valence 
band. The electron and hole are oppositely charged and 
may form a bound state, the exciton, with the absence or 
presence of a ground state exciton within a dot forming 
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our qubit basis (|0) and |1) respectively). For excitons, 
we must consider an electron-hole pair Hamiltonian 



H = H e + Hh 



47re(r e - r h )|r e - r h 



+ E gap , (7) 



where H e and Hh are given by Eq. ^ with the appro- 
priate effective masses and potentials; E 9ap is the semi- 
conductor band-gap energy, and e(r e — rh) is the back- 
ground dielectric constant of the semiconductor. Wc 
shall consider the simplest case of e(r e — rh) = eoe r , i-e. 
the relative permittivity e r is independent of (r e — rh). 
The intra-dot energy shift due to the Coulomb term 
H e h = e 2 /47reo e r| r e — r h| is a small contribution to the 
total energy and we treat it as a first-order perturba- 
tion. This strong confinement regime treats the electron 
and hole as independent particles with energy states pri- 
marily determined by their respective confinement po- 
tentials i*» It is valid for small dots with sizes less than 
the corresponding bulk exciton radius ao (~ 35 nm for 
InAs, ~ 13 nm for GaAs). For more sophisticated treat- 
ments of the calculation of excitonic states see, for exam- 
ple, Refs. fdirect-diagonalization) . l2ll (psuedopotential 
calculations), and I25I (variational methods). In Ref. I21I 
it was found that a simple perturbation method was in 
good agreement with a self-consistent field approach. 

We construct an antisymmetric wavefunction repre- 
senting a single exciton state given by 



tfi = A[V4(ri,(7i),Vm(ra,a-a)] 



(8) 



where r and a are position (from the centre of the 
dot) and spin variables respectively, n and m label the 
quantum states, and A denotes overall antisymmetry. 
Here, one electron ^(ri,<7i) has been promoted from 
the valence band into a conduction band state whilst 
V-Vn( r 2, 02) represents a state in the valence band. Tak- 
ing the Coulomb matrix element (^F\H e h\^i) between 
the initial state ^1 above and an identical state \I/f (in 
effect coupling an electon and hole via the Coulomb op- 
erator) leads to two terms^SiSl the direct term 
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and the exchange term 



Exch 
IF 



= ±- 



tp'*(ri)tp m (Tx) 



x V4( r 2)VC( r 2)dridr 2 . 
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(10) 



The sign of the exchange term is determined by the 
symmetry of the spin state of the two particles; with the 
perturbation H e h being positive, triplet spin states give 
negative exchange elements whereas singlet spin states 
give positive values. We shall now show how to calculate 
the direct electron-hole Coulomb matrix element on a 
single dot where n and m are both taken as ground states. 



The exchange interaction is much smalle»2i*2& and we 
shall not consider it here. 

If we consider identical potentials in all three direc- 
tions then we can use the spherical symmetry to derive 
an analytical expression for the direct Coulomb matrix 
element, which we call M e h- For d x = d y = d z = d, Eq. [21 
may be written in spherical polar coordinates as 



(r) = 



1 



3/2 



Substituting into Eq. El leads to 
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M eh = 



x exp 
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fl2) 



where the contribution of the Bloch functions U(r) has 
been neglected^ We now express l/|ri — r 2 | in terms of 
Legendre polynomials as2& 



n -r 2 



£E,^o(£)'fl(«»< 



for n > T2 
for r\ < r 2 . 



(13) 

Substituting this into Eq. 1121 and integrating over polar 
angles leads to 
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(14) 



where use has been made of the orthogonality relations 
of Legendre polynomials. The integrations now give us 
the following expression for M e h'- 
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In a similar manner, we may also approximate the be- 
haviour of M e h in the presence of an external electric 
field. For a constant field applied to the dot, the poten- 
tial in the field direction (for simplicity say z, although 
the spherical symmetry we assume means all three direc- 
tions are equivalent) becomes 



V(z) 1 — * V(z) + qFz, 



(16) 



where q = — e for conduction band electrons, q = +e for 
holes, and F is the electric field strength. Substituting 
this into the Schrodinger equation for the ^-component 
leads us to a new Schrodinger equation that has the same 
parabolic potential form: 
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{z')=E'ct>{z') (17) 
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with 



z' e = z- eF/c C:Z 
z' h = z + eF/c h , z 
E' i = E+ (eF) 2 /2c v 



for electrons, 



for holes, 



(18) 



e, h. 



Therefore, electrons and holes are displaced in opposite 
directions and their envelope functions are the same as 
Eq.[|J]with z replaced by z'. The simplicity of the change 
in envelope function with applied electric field is a great 
advantage of the parabolic well model, although it should 
be pointed out that this same simplicity implies that the 
charges can continue separating indefinitely with applied 
field strength and is therefore unrealistic at very high 
fields. 

Again, in spherical polar coordinates, the envelope 
functions in the presence of a field may be written as 



<Mr) 
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for electrons, and 
1 
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exp 
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(19) 



(20) 



in a up to the term in a 3 and integrate over r\. Keep- 
ing only the terms up to (d e a) 2 order in the resultant 
expressions gives us an estimate for the suppression of 
the electron-hole binding energy as an external field is 
applied: 
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(25) 



which reduces to Eq. ^|at F = 0. 

Equations 1181 and 1251 imply a quadratic dependence 
of the Stark shift (change in exciton energy) on the ap- 
plied electric field. This has been observed experimen- 
tally in a range of QD systems including InGaAs/GaAs, 29 
GaAs/GaAlAsr2i and CdSe/ZnSe^i Furthermore, a the- 
oretical study of an eight-band strain dependent k • p 
Hamiltonian has shown that the quadratic dependence 
of the ground-state energy on applied field is a good ap- 
proximation for largely truncated self-assembled quan- 
tum dots, although the approximation becomes worse as 
the dot size increases in the growth direction. 32 



for holes, where k is the unit vector in the ^-direction. 
This time, substituting into Eq. pleads to 
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We proceed as before, again making use of Legendre poly- 
nomials and their orthogonality relations, and integrate 
over r 2 to get 



M eh = ^ / exp 



erf 



ri 



a J \ d-U \dh 
x [exp (ria) — exp (— r\a)] dri 



(24) 



For a <§; l/d e (valid up to fields of order 10 7 V/m for the 
small dots considered here), we expand the exponentials 



III. A SIGNATURE OF FORSTER COUPLED 
QUANTUM DOTS 

A. Hamiltonian 



We have now characterized the single particle electron 
and hole states within a simple QD model, as well as 
accounting for the binding energy due to electron-hole 
coupling within a dot when estimating the ground state 
exciton energy. In this section we shall consider exci- 
tons in two coupled QD's and the Coulomb interactions 
between them. More specifically, we shall derive an ana- 
lytical expression for the strength of the inter-dot Forster 
coupling. We shall show that this coupling is, under cer- 
tain conditions, of dipole-dipole typ o 3 i 10 and that it is 
responsible for resonant exciton exchange between ad- 
jacent QD's. This is a transfer of energy only, not a 
tunnelling effect. We are concerned in this paper with 
bringing excitons within adjacent QD's into resonance. 
As the Appendix shows, single particle tunnelling is only 
significant when the energies of the states before and af- 
ter the tunnelling event are separated by less than the 
tunnelling energy. This is a different resonant condition 
to the one considered here and is not fulfilled by the dots 
over the parameter ranges explored. 

Following Ref. |^ we write the Hamiltonian of 
two interacting QD's in the computational basis 



5 




FIG. 2: Schematic diagram of the interacting two-dot system. 



{|00>, |01>, |10), |11>} as (h=l) 



H = 














UJq + LU 2 


V F 





v F 


w + 
















ujq + uj± + uj 2 + Vxx ) 



(26) 

where the off-diagonal Forster interaction is given by Vf , 
and the direct Coulomb binding energy between the two 
excitons, one on each dot, is on the diagonal and given 
by Vxx^ The ground state energy is denoted by ojq, and 
Auj = u>i — ll>2 is the difference between the excitation 
energy for dot I and that for dot II. These excitation 
energies and inter-dot interactions are all functions of 
the applied field F. The energies and eigenstates of this 
four-level system are given by 



-E-00 — ^0, 

= LJ ■+ 

E + = uj 4 
En = w - 



LUl - 
Ll>1 



A), 

4p(l - A), 
w 2 + Vxx, 



|*oo) 
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l*u) 



= |00) 

= ai|10)-a 2 |01) 
= ai|01)+a 2 |10) 
= 111), 



where A = y/1 + 4(V f /Acj) 2 , a\ 
a 2 = sgn(VFAw)>/(A 



= vV^ 

1)/2A for | Aw | > 



(27) 
and 



1)/2A 
0. We can 

see that Vf may cause a mixing of the states |01) and 
1 10) with the result that |'F_) and |^+) can now be en- 
tangled states. It is also straightforward to see that an 
off-diagonal Forster coupling does indeed correspond to a 
resonant transfer of energy; if we begin in the state 1 10} 
(exciton on dot I, no exciton on dot II) this will natu- 
rally evolve to a state 1 01} (no exciton on dot I, exciton 
on dot II), in a time given by 7r/(2Vp), through the maxi- 
mally entangled state 2 -1 / 2 ( 1 10) + i|01)). An analogous 
behaviour is expected for the initial state |01). 



B. Analytical model of the Forster interaction 

We shall now calculate the magnitude of the off- 
diagonal matrix elements in Eq. 1261 for the parabolic po- 
tential model. We shall see that the behaviour of the 
Forster interaction due to changes in dot size, composi- 
tion, separation, and applied electric fields may be pre- 
dicted by such an analytical model. We begin by calcu- 



lating the form of the matrix element in the dipole-dipole 
approximation. 

The matrix element we require is that of the Coulomb 
operator between two single exciton wavcf unctions, one 
located on each of the two dots. We take our initial state 
as representing a conduction band state in dot I, and a 
valence band state in dot II: 



*I = - 4 [V , n( r Xj°'l))V ! m(r2,0'2)] 



(28) 



For our final state, we must have a valence band state in 
dot I and a conduction band state in dot II, given by: 



* F = A[^ n {v 1 ,a 1 ),il)' m {Y 2l a 2 )} ■ 



(29) 



The positions ri and r 2 are now defined from the centres 
of dot I and dot II respectively (see Fig. [2J , and not from 
the same point as in Eq. [HI where only a single dot was 
considered. Note that the time ordering of \&i and 
is irrelevant since, as described in Sec. II I II the resonant 
energy transfer process is reversible. 

Therefore, the direct Coulomb matrix element between 
these two states gives us 



Vf = 



C(ri)^(ri) 



1 



47reo£r 
x VCMVCMdrxdra, 



iR + ri-ral 



(30) 



where we have explicitly included the inter-dot separation 
R. For |R| ^> |ri — r 2 |, which is valid as long as the 
characteristic sizes of the wavefunctions, dj, are small 
in comparison to |R|, we can follow the procedure of 
Dexte»i^ and expand the Coulomb operator in powers of 
( r i/ 2 /!R) up to second order. Taking the matrix element 
between and ^f leads to 



V F 



47re e r i? 3 



(ri) • (rn) 



R 2 



«rr).R)«rn).R) 



(31) 

We assume the dots are sufficiently separated for there 
to be no overlap of envelope functions between dot I and 
dot II. Therefore, we do not consider the exchange term 
of this Coulomb interaction. The integrals 



(r I )=/<*(r 1 )r 1 ^„(ri)dr 1 , 
(r n ) = / ^(r 2 )r 2 ^(r 2 )dr 2 , 



(32) 



are taken between an electron and hole ground state cen- 
tered on dot I and dot II respectively. Remembering that 
our wavefunctions are a product of an envelope function 
4>(r) and a Bloch function U (r) we can make use of their 
different periodicities to write2& 
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47reoe r i? 3 
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(33) 
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FIG. 3: Dependence of the overlap integral of a single 
parabolic QD as a function of the in-plane (x- direction) elec- 
tric field strength. Suppression of the Forster interaction re- 
sults from the reduction in electron-hole overlap for increasing 
fields. A 162.5 meV potential at x = 3 nm from the dot centre 
for both electrons and holes is used, giving c = 0.00579 J/m 2 
(see Fig. . 



The overlap integrals are defined as 



O 



(r)</> h (r)dr, 



with Oi/n referring to the overlap of the envelope func- 
tions for dot I or dot II respectively (each having a max- 
imum value of unity), and the inter-band dipole matrix 
elements are defined as 



U e (r)rU h (r)dr, 



(35) 



with d cv (i/n) referring to dot I or dot II respectively. 

We shall not calculate the values of d cv here as they 
are commonly measured experimental quantities (see also 
Ref. |2(| for a simple model) and, once the dot materi- 
als have been chosen, are constant contributions to the 
Forster interaction strength. However, the calculation of 
0\/u is vital in determining the effects of dot size, shape, 
and applied electric fields on the strength of the inter-dot 
interaction. 

We take the parabolic solutions in Cartesian coordi- 
nates from Eq. and also include the effect of a lateral 
electric field, which is important in determining how to 
suppress the interaction when required, or bring two non- 
identical dots into resonance. As before, we shall assume 
that the electric field affects only the envelope function 
part of the wavefunction; this is valid in the regime where 
the electric field never becomes so large that the envelope 
function varies on the unit cell scale. This is equivalent 
to saying that the envelope functions can be decomposed 



into a superposition of crystal momentum (k) eigenstates 
near the band edges, where the Bloch functions are ap- 
proximately independent of k. 

For a constant field in the lateral direction (say, x) the 
overlap integrals Oi/n are straightforward to calculate. 
Again, for clarity, we take wells identical in all three di- 
rections for both electrons and holes {d x = d y = d z = d), 
leading to 



Mr) 
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x exp 



3/2 



exp 



(x - eF/c e f 
2dl 



(y 2 + z 2 ) 

2di 



(36) 



for electrons, and 
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( 1 



x exp 
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exp 
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(y 2 + z 2 ) 



2d l 



(37) 



for holes. Substituting into Eg. 1341 and integrating results 
in 



O = 



2d e dh 



dl 



dl 



3/2 



exp 



e 2 F*(c e + c h f 



2cH(dl 



dl) 



(38) 



Therefore, in zero applied field the overlap depends only 
(34) on the ratio (2d e dh/(dg + d%j) S . It is worth noting that 
if we had chosen an infinite square well potential in the 
growth (z) direction and parabolic wells in the x- and 
^-directions (as is common in the literatureii) then the 
zero field value would be (2d e dh/{d% + d\)) with the field 
dependence being exactly the same as in Eq. 1381 

Figure shows the suppression of the overlap integrals 
by an in-plane electric field (^-direction) and therefore 
the suppression of the Forster interaction itself. How- 
ever, as we shall see in the next section, this does not 
rule out its observation in coupled dot systems that are 
tuned to resonance with an external applied field. Indeed, 
if the inter-dot interaction matrix element is relatively 
large in zero applied field, then interesting anticrossing 
behaviour should be observed in the energy spectrum as 
the system is tuned through resonance. We also note 
that a suppressed Forster coupling can be of benefit to 
the exciton-exciton dipole interaction quantum compu- 
tation schemes^ as it ensures an almost purely diagonal 
interaction between adjacent QD's. 

Taking measured values for the transition dipole mo- 
ment e(r) allows us to estimate the magnitude of Vp be- 
tween two stacked dots. In CdSe QD's a value for e(r) 
of up to 5.2 eA has been reported^ while for both In- 
GaAs/GaAs and InAs/InGaAs QD's values of approx- 
imately 5-7 eA have been measuredi^24 Considering 
Eq. |^ with a value of 6 eA for e(r), e r — 12 (for In- 
GaAs/GaAs), and inter-dot spacing R = 5 nm, we ob- 
tain an estimate of 0.69 meV for the Forster coupling 
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FIG. 4: Forster interaction strength as a function of dot sep- 
aration for two identical dots, with e r — 12. Three different 
values of e{r) are shown: 7 eA (solid line), 6 eA (dashed line), 
and 5 eA (dotted line) 



energy Vf, certainly large enough to be observed experi- 
mentally. This corresponds to a resonant energy transfer 
time of picosecond order and is therefore interesting as a 
coupling mechanism for performing quantum logic gates, 
as it is well within the nanosecond dephasing times^22i2£ 
expected for excitons within QD's (see Sec.E]for further 
discussion). In Fig.0]the 1/R 3 dependence of the Forster 
interaction strength is shown for various values of e(r). 

In the next section we shall discuss a signature of 
the Forster interaction that would be observable through 
photoluminescence measurements. 



the resonance, an anticrossing of the energy levels should 
be observable through photoluminescence measurements. 
However, the transition from the antisymmetric state to 
the ground state is not dipole allowed on resonance and 
should also display a characteristic loss of intensity close 
to the resonant condition (see Sec. II VJl . 

From Eq. ^| we see that an external field reduces the 
energy E for both electrons and holes, with the shift be- 
ing greater for bigger dots. We therefore consider two 
coupled dots of different material concentrations, with 
one of slightly greater dimensions and having a larger 
band-gap, and in Fig. |5£a) show that the single dot en- 
ergy levels cross (for our choice of parameters) as an elec- 
tric field is applied. Such a situation is plausible for sys- 
tems such as InGaAs where dot layers of varying indium 
content, and hence varying band-gap, may be grown2& 
and it applies directly to the parameters chosen here 
(see also Sec. IVI below). Diagonalising the Hamiltonian 
(Eq. I2rj|) as the held strength F varies gives us a model 
prediction for the behaviour of the energy levels £L and 
E+ as shown in Fig. |5[b) , where an anticrossing is ob- 
served at a field of approximately 7.5 x 10 6 V/m. Since 
u>x, u>2, and Vf are all functions of F (as is also shown 
in Eqs. El an d I38[) there is large scope for finding 
parameter regimes with interesting behaviour. 

An analytical expression for the field strength at reso- 
nance can be calculated from the condition Alu = 0. The 
total energy of each dot is given by 

E = E*» + f ( JT + JT) - £2! - M eh> (39) 

2 VV m e V m hJ C 

for identical potential wells in all directions for both elec- 
trons and holes (c x = c y = c z = c). When dot I and dot 
II are resonant 



Ei - En = 



C. Anticrossings: A signature of Forster coupling 

If we consider again Eq. |23we can see that and 
\^ +) have a range of forms depending upon the values of 
a\ and ai. For example, if Aw 3> Vf, then A ~ 1 which 
leads to |ai| ~ and |oa| — 1- Therefore the states \^~) 
and |\& + ) are given by 1 01) and 1 10) respectively and there 
is no mixing of the computational basis states. The only 
way to couple two dots in this case is via the diagonal 
interaction Vxx- However, for two dots coupled by Vf 
at resonance (Aui = 0) we can see from Eq. |23 and 
by using A = ^{Aoj 2 + 4Vp 2 )/Aw 2 , that |a x | = |a 2 | = 
1/ y/2, i?iowcr = wo + ui — \ Vf\, and E highcr = uj + w\ + 
\V ¥ \. Furthermore, the two eigenstates 2 _1 / 2 (|01) + 1 10) ) 
and 2 _1 / 2 (|10) — 1 01) ) are both maximally entangled and 
separated in energy by 2Vp- 

Interestingly, we should be able to move between these 
two cases by bringing two initially non-resonant coupled 
dots into resonance, for example by the application of a 
static external electric field. By taking the dots through 



{E 9 l ap - E^P) + — 



1 












i 




ClI 



and therefore 



F 




(E( ap - Eg*) 



'ml 



2tt 3 / 2 



e e r 



«+<) 1/2 «+< : ) 1/2 



(41) 
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FIG. 5: (a) Single dot energies as a function of applied electric field. Parameters: ml — 0.04mo, m* h = 0.45rao, e r = 12, and 
ci = 0.00579 J/m 2 (corresponding to a potential of 162.5 meV at a distance of 3 nm from the dot centre) for dot I, for both 
electrons and holes, and cn = 0.03414 J/m 2 (corresponding to a potential of 515.7 meV at a distance of 2.2 nm from the dot 
centre) for dot II, for both electrons and holes (see Fig. 0. E 9ap is taken as 1.1 eV for dot I and 0.394 eV for dot II. (b) 
Energies E- and E+ of the coupled dot system demonstrating anticrossing at a field of approximately 7.5 x 10 6 V/m. Here, 
Vf has a magnitude of 0.20 meV at zero field, with d cv = 7 eA and R = 7 nm. 



where 



2e 2 



1 



1 



1 



37T 3 / 2 £o e r {d 

2e 2 



1 



<) 3/2 ci 
1 



(42) 



and the states |\&_) and |\&+) should be maximally en- 
tangled at this value of F (= 7.547 x 10 6 V/m with the 
same parameters as for Fig. |SJ), with an energy separa- 
tion equal to 2Vp (at this field) as stated earlier. Clearly, 
the experimental observation of an anticrossing as shown 
in Fig. Elb) would be an extremely strong indication of 
Forster coupling between two dots, and also a first indi- 
cation that entangled states are being produced. 



IV. DECAY RATES AND ABSORPTION 

We have seen in the previous section that an anticross- 
ing in the energy level structure of two coupled QD's 



provides a signature of the Forster interaction, which 
should be observable through photoluminescence mea- 
surements. However, the scenario considered thus far 
is idealized in that there is no coupling of the two-dot 
system to the external environment. Any emission (or 
absorption) lines associated with coupled-dot transitions 
will also be broadened due to emission ^ and scattering 
and pure dephasing processes due to exciton-phonon in- 
teractions^ Experimentally, Bayer and Forchel 36 have 
shown that decay processes are dominant at very low 
temperatures (~ 2 K), while other studies have demon- 
strated good Lorentzian fits to the photoluminescence 
lineshapes at 5 Hence, we shall limit the discus- 

sion here to spontaneous emission (decay) processes. 

The spontaneous emission rate for a two-level system 
(single dot) interacting with a single radiation mode is 
usually calculated by considering a quantum mechani- 
cal description of the radiation field (see, for example, 
Ref. l4l|) , although approaches which consider a classical 
light field also exists For a dot surrounded by material 
of approximately the same relative permittivity e r the 
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decay rate becomes^ 



Od c 



3irc 3 heo 



(43) 



where huio is the energy difference of the two levels un- 
der consideration. For a typical InGaAs dot we take the 
parameters uio = 1.3 eV, Od cv = 6 eA, and e r = 12 to 
give T sp = 1.04 x 10 9 s _1 or a decay time of T^ecay = 
l/r sp = 964 ps, comparable with experimentally mea- 
sured exciton lifetimes in this systemi2£iSii2£ 

Here, we are primarily interested in the properties of 
two interacting dots which form the four-level system 
considered in Sec. IIII Al The various decay rates be- 
tween each level may be calculated in the same manner 
as for the two-level system previously considered, provid- 
ing that the changes in transition dipolc moments due to 
the interaction are properly accounted for. We will then 
be able to predict the typical linewidths that would be 
observed in experimental measurements of these transi- 
tions. We characterize the dipole operator in the com- 
putational basis according to which dot the transition 
occurs within: 




5 10 

Field (MV/m) 
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FIG. 6: Spontaneous emission rates of the coupled-dot energy 
levels of Fig. |^1 The dashed curve corresponds to the upper 
curve in Fig. (b) ; the solid curve corresponds to the lower 
curve in Fig. [51(b). 
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Transitions such as 1 1 1 ) — > 1 00) have zero dipole moment 
since the corresponding integral is zero due to the orthog- 
onality of valence and conduction band wavefunctions on 
each dot: 

(r) = J <(r 1 X(r 2 )(r 1 +r 2 )<(r 1 )^(r 2 )dr 1 dr 2 
= / K( r i) r iV4( r i)dri / Vm( r 2)^m( r 2)dr 2 



+ / C(ri)^(ri)d ri J ^(r 2 )r 2 ^(r 2 )dr 2 



0. 



(45) 



As a result of Eg . 1441 we may express the dipole moments 
for general transitions such as <z|01) ± 6 1 10) — > |00) by 

(r) 0l6 = o<00|r|01> ± 6<00|r]10) = aO n d cv(I1) ±bO I d cv(I) , 

which may then be inserted directly into Eq. 1431 along 
with the correct frequencies, to give the corresponding 
decay rates. In Fig. [5] we plot T sp for the two energy 
curves of Fig. |3] (b) from Eqs. E7J ESI an d ED an d with 
the same parameters as Fig. [S] A special case occurs for 
identical dots at the anticrossing. Here, the symmetric 
eigenstate 2 _1 / 2 (|01) + 1 10)) has a transition dipole mo- 
ment to the ground state of \/20d cv (0\ = On = O, 
d C v(i) = d cv (n) = d cv ) and hence a decay rate of twice 
that expected for a single dot. However, the antisymmet- 
ric eigenstate 2 _1 / 2 (|01) — 1 10) ) has no transition dipole 



moment and consequently no spontaneous emission rate 
in the dipole approximation. This is otherwise known as 
a "dark" state and any spectral line corresponding to this 
transition will display a characteristic narrowing and loss 
of intensity as the anticrossing is approached. 

This effect may be studied in more detail by consid- 
ering the absorption lineshape of each transition (within 
the rotating- wave approximation)!^^ 



(T s p/2) 



c^heo (wio - + (r*P/2) 2 



(47) 



where r\ is the refractive index. Here, the only line broad- 
ening which is accounted for is due to spontaneous emis- 
sion. This leads to absorption lines with a Lorentzian 
dependence on frequency, and a full width at half maxi- 
mum given by T sp . Although other mechanisms may also 
broaden the lines, for example "pure" dephasing due to 
exciton-phonon interactions as mentioned earlier, these 
processes can usually be reduced, in our case by cooling 
the system^ However, it is difficult to reduce the spon- 
taneous emission rate of a given transition. Therefore, 
the linewidth T sp is the minimum achievable from any 
standard dot sample and hence it is vitally important to 
ensure that any effects we wish to observe will not be 
masked by its presence. 

Plotted in Fig. [7J are the absorption spectra of the 
two energy levels in Fig. GJb) (E + and E-) at fields of 
0, 5, 7, 10, 15 MV/m respectively, calculated from Ed. 1371 
with the same parameters as for Fig.[S]and -q = 3.46. The 
spontaneous emission rates are calculated from Eq. EHfl 
and the peaks have been artificially broadened by a fac- 
tor of 10 to exaggerate their characteristic features in 
the changing applied field. As the field is increased the 
two peaks shift to lower energies with the initially higher 
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FIG. 7: Series of simulated absorption spectra of the energy 
levels in Fig.^b) at fields of F = 0, 5, 7, 10, 15 MV/m. These 
lines have been artificially broadened by a factor of 10. 



energy line (corresponding to E + ) shifting by a greater 
amount so that their separation reduces. The width of 
the lower energy line (£L) increases as the anticrossing 
point (F = 7.5 MV/m) is approached signifying its in- 
creasing decay rate (see Eq. 14611 . On the other hand, 
the width of the E+ line decreases as it approaches a 
dark state, and the area underneath the curve is re- 
duced. This would correspond to a lowering of intensity 
of this line in a photoluminescence experiment. We can 
see that the separation of the lines due to the Forster 
interaction (2Vp ~ 0.4 meV) is well resolved in the pres- 
ence of radiative broadening at low temperatures, and 
should be so even if the lines have extra broadening due 
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to exciton-phonon interactions. In fact, sharp emission 
lines of approximately 0.1 meV width have been obtained 
from single InGaAs QD's at a temperature of 100 K, 36 
indicating that up to this temperature at least, the an- 
ticrossing effect should still be observable. At very high 
fields (F = 10 — 15 MV/m) beyond the anticrossing point 
the two lines once again become well separated and even- 
tually have similar widths, indicating that the states 1 01) 
and 1 10) are now only weakly coupled. 



V. CONNECTION TO QUANTUM 
INFORMATION PROCESSING 

The experimental observation through photolumines- 
cence of an anticrossing of the type above would be a 
significant step towards proof-of-principle experiments; 
although it could be a difficult experiment to perform, 
this method may well yield results more quickly than an 
attempt at coherent control on dot systems coupled in 
this way. 

The main question here is the feasibility of bringing 
two dots into resonance using a static external electric 
field. As has been mentioned above, and can be seen in 
Fig. OJa), a larger dot experiences a larger shift in its 
energy levels due to the applied field than a smaller one. 
However, all other parameters being equal, a larger dot 
also has slightly lower energy at zero field (which is not 
the case in Fig. [5J. Therefore, a way is needed of in- 
creasing the initial energy of the larger dot relative to 
the smaller one. This could be realised by using layers 
of different materials (or material concentrations) to al- 
ter the band-gap within each dot; other methods such 
as exploiting different dot geometries or applying a local 
strain or electric field gradient should also be explored. 
In Rcf. 29 field gradients close to 20 (MV/m)//xm were 
generated, with Stark shifts of approximately 2 meV ob- 
tained in a field of 0.2 MV/m. Hence, similar dots of 
1 — 2 meV initial energy separation, and placed 7 run 
apart, could be brought into resonance by this method. 
Nitride QD's could also offer a promising approach since 
their strong piezoelectric fields allow the possibility of 
an external field shifting their energy levels towards each 
other. 

We have also shown that it should be possible to engi- 
neer nanostructures such that the off-diagonal Forster 
interaction between a pair of QD's is of the required 
strength to make it interesting for quantum computation. 
Once a measurement of this coupling strength is made, 
the next logical step is to attempt to controllably entan- 
gle the excitonic states of two interacting dots, leading on 
to a demonstration of a simple quantum logic gate such 
as the controlled-NOT (CNOT). Although, on resonance, 
the states 1 01) and 1 10) naturally evolve into maximally 
entangled states after a time 7r/(4Vp), their initializa- 
tion requires the inter-dot interaction to be suppressed. 
Furthermore, the generation of a logic gate such as the 
CNOT requires single qubit operations on both dots, as 
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well as periods of interaction. 

By switching to a pseudo-spin description of our ex- 
citonic qubit we can immediately consider a previously 
known operation sequence for the realization of a CNOT 
gate. Defining | | z ) = |0) and | [ z ) = |1), we can see 
from Eq.[23that the off-diagonal terms can be expressed 



Hp = — (a Xl a X2 + <J yi (jy 2 ) 



with 



1 

1 



and er„ 



-i 

1 



(48) 



(49) 



being two of the Pauli spin matrices. This XY type 
Hamiltonian has been studied in the literature for various 
systems^^i and if the two interacting qubits are left for 
a time t = tt/(2Vf) solely under its influence then an 
iSWAP gate will be executed^ 



iSWAP 



n o o o\ 

i 

i 00 

\0 1/ 



(50) 



in the basis {|00), |01), |10), |11)}. Two iSWAP opera- 
tions may be concatenated with single qubit operations 
to form the more familiar CNOT gate: 



CNOT = ( - 



(iSWAP)® 



(51) 



where (±7r/2); m are single pseudo-spin rotations of ±7r/2 
about the / axis of spin m, for I — x,z and m = 1,2. 
Schuch and Siewerl^ have also shown that the CNOT 
and SWAP operations may be combined when using an 
XY interaction to produce more efficient quantum cir- 
cuits. Furthermore, the iSWAP operation is an entan- 
gling gate and is therefore sufficient for universal quan- 
tum computation provided that fast local unitary opera- 
tions are available. In fact, for systems exhibiting an XY 
interaction, the iSWAP operation constitutes the natu- 
ral gate choice when implementing efficient quantum cir- 
cuits. 

To perform a gate such as the CNOT outlined above 
we must be able to control the interaction between our 
two qubits so that we can effectively switch it off for the 
duration of the single qubit manipulations. For the case 
of excitonic qubits, coupled via the Forster mechanism, 
the most sensible way to proceed is to consider two ini- 
tially non-resonant QD's with negligible energy transfer. 
Single qubit operations can then be achieved with ex- 
ternal laser pulses by inducing Rabi oscillations within 
each dot 4^ As each dot will have a different excitation 
energy, we may address them individually by choosing 
the appropriate frequency. Two periods of free evolu- 
tion under the interaction Hamiltonian (Eq. I48|) are also 



required; applying a suitably selected detuned pulse to 
both non-resonant dots will bring them into resonance 
via the optical Stark effecti^iLSi We then allow resonant 
energy transfer to occur for a time t — ir/ (2Vf) producing 
an iSWAP operation. The detuned pulse is then stopped 
and single qubit manipulations may be induced as before. 

Figure EJb) provides a nice visualisation of the whole 
process. We must non-adiabatically switch between the 
two regimes of zero field, where the dots are effectively 
uncoupled, and the resonant point where the dots inter- 
act. It is our hope that this is achievable through the 
optical Stark effect, and we speculate that this all opti- 
cal approach may have the potential to allow gates to be 
performed well within the limits set by the nanosecond 
dephasing times experimentally observed. 

To summarise, we have analytically calculated the 
magnitude of the Forster energy transfer between a pair 
of generic QD's and investigated its effect on their energy 
level structure. We have proposed a simple experiment 
which provides a signature of the interaction and an esti- 
mate of its strength, and have also discussed its possible 
application to quantum information processing. 
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APPENDIX: TUNNELING 

To outline the effect of electron and hole tunnel- 
ing on the exciton states used in this paper, we con- 
sider here the Hamiltonian for an electron-hole pair 
in a double-dot. The basis we use is constructed of 
products of the electron and hole single-particle states 
{|eifti), \eihu), \enhi), \euh u )}, which gives 



H = 



(E ei hi th t e Vf 
th E ei h u t e 

t e ^en hi th 
\ Vf t e th Sen hi 



(A.l) 



where E en h m = E Cn + Eh m — M en h m , with n, m = I, II 
for dot I and dot II respectively. M &ri h m is the direct 
Coulomb binding energy between the electron and hole 
on dot n and m respectively, and the band-gap energy has 
been absorbed into the electron energy E Cn by setting the 
energy zero to be at the top of the valence band. Vf is 
the Forster interaction strength, and i e (/i) is the electron 
(hole) tunneling matrix element. 
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We consider first the simple case of two identical dots 
coupled to one another; this means setting E ei — E en = 
E e , Ehj — Eh n = Eh, M ei /ti = M eu i lu = M e h, and 
M ei h u = Me u hi = M' eh . Su btracting E e + E h - M eh 
from the diagonal of Eq. IA.1I gives 



H = 



■ 


th 


t e 


Vf 


th Meh 


-K h 





te 


t e 


M eh 


-M^h 


th 


Vf 


te 


th 






(A.2) 



We would like to isolate the {|ei/ii), |en/in)} subspace, 
as it is composed of single exciton states on each of the 
two dots. These states are exactly the ones that are rel- 
evant for the computational basis introduced in Sec. [I] 
Any leakage from this subspace, potentially due to tun- 
nel couplings to the states \e\h\i) and \e\\h\), could be a 
source of error for the signature and schemes presented 
here and in Refs. l8l26l and must be minimized. However, 
under the condition 



(A.3) 



\M eh -M' eh \ » \t e \,\t h \, 



we may use degenerate perturbation theory on Eq. I A. 21 
to give 



ti+ti 



M, 



-Kh 
2t.t h 



2t e t h 



it+i 



-M' 



(A.4) 



M, 



-Ml,. 



M e , 



-M' 



in the {\e\h\), \euh\i)} subspace. Hence, the states \e\h\) 
and \e\ihii) are still resonantly coupled in the presence 
of tunneling as long as Eq. I A. 31 is satisfied. These 
conditions are better satisfied as the inter-dot separa- 
tion increases (tunneling elements consequently reduce, 
as does M' eh so that \M e h — M' eh \ becomes larger), 
and as dot confinement increases (tunneling elements 
reduce, M e h increases so that \M e h — M' eh \ again be- 
comes larger). Furthermore, corrections to the eigen- 
states \x±) = 2 _1 / 2 (|ei/ii) ± \enhu)) due to mixing with 



states outside the subspace will be small since they are 
weighted by factors of t e ( h )/{M' eh — M e h), to first order, 
from the perturbation theory. 

The regime in which Eq. I A. 21 is valid is not necessar- 
ily the ideal one for minimizing the effect of tunneling, 
while exploiting resonant exciton interactions, as two 
non-identical dots may also be brought into resonance 
(see Section fill C|l . In this case, E ei + E^ — M ei hi — 
E en + Eh u — M eu h u = E on resonance. Subtracting E 
from the diagonal of Eq. IA.1I gives 



H 



where AE q 

M eihl - M eih 






t h 


t h AE h 
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te 
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te 


= E in 
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AE e + AM e 
t h 



(A.5) 



£L; T , for 



AM, 



M 



eihi 



M eu hi- 



e,h, and AM h = 
The dots must 



now satisfy the modified condition 

min(\AE h + AM h \, \AE e + AM e \) > \t e \, \t h \, (A.6) 
in order for tunneling to be neglected, with the unwanted 
states | ei/z-n) and \euhi) weighted by a factors of magni- 
tude 



e(h) 



AE h + AMh 



(A.7) 



and 



t 



e(h) 



\AE P + AM P 



(A. 



to first order in a perturbation expansion. Again, tunnel- 
ing will be suppressed as dot separation and confinement 
increases. 
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